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1 Introduction 

In P a mathematical model of an anomalous diffusion in a space was suggested. This model 
origins in an investigation of processes in complex systems with variable structure: glasses, 
liquid crystals, biopolymers, proteins and etc. 

In this model a coordinate of a particle has stable distribution (not normal one). As a 
result a density of its distribution function satisfies an analog of diffusion equation in which 
second derivative by coordinate is replaced by partial derivative. 

In this paper the anomalous diffusion with periodic initial conditions on an interval with 
reflecting edges is considered. Such problem is important for example in technical mechanics 
for an analysis of fuel mixing in straight flow engine [2] too. As A. A. Borovkov suggested [B] 
special probability methods are developed to analyze such a diffusion. 

Main results of the paper may be represented as follows. Suppose that T n (a,r,k) is a 
characteristic time of the anomalous diffusion with a parameter a, < a < 2 on fc-dimension 
interval [— r, r] h with n-periodic (by each coordinate) initial conditions then 

T 1 (a,r,l)=T 1 (a,l,l)r a , ^r^^ = r a ~ 2 , T n (a, 1, 1) = H(a, 1, l)/n«\ 

Ti[a, 1, l)Ti(2,r, 1) 

So the anomalous diffusion on the interval [— r, r] for r > 1 works faster and for r < 1 - slower 
than normal one. If k > 1 then it is possible to choose n-periodical (by each coordinate) 
initial conditions so that for the normal diffusion 

T n (2,l,k)=T 1 (2,l,l)/kn 2 . 



2 Model of Anomalous Diffusion on Straight Line 

Suppose that y(t), t > is homogenous random process with independent increments 
and an initial condition y(0) = 0. The difference y(t) — y(r), t > r > has symmetric 
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stable distribution on the straight line (—00, 00) with a parameter a, < a < 2 and a 
characteristic function 

M exj>(iu[y(t) - y{r))) = exp(-(t - r)\u\ a ). 

In P it was shown that distribution density p t = Pt{u) of the process y(t) in a moment t > 
is the generalized solution of the following differential equation with fractional derivatives 

Here d a p t (y) / dy a is a— fractional derivative of the function pt{y), < a < 2; <5(i) are 

delta-functions by variables y, t accordingly 

3 Construction of Anomalous Diffusion Model 
on the Interval [—1, 1]. 

Each realization of random process y{t), t > may be considered as a curve T on the plane 
(y, t) and may be represented in parametrical form: y = yij), t = t{r) = r, r > 0. Suppose 
that the curve T is reflected from the lines y = 1, y = — 1. For this aim represent the plane 
(y, t) as a transparent and infinitely thin sheet of paper with the curve T. Bend this sheet 
of the paper along the lines y = ±1, y = ±3, . . . into transparent strip — 1 < y < 1 with 
fragments of initial curve T. As a result T is transformed into the curve 7 : y = Y(t), t = 
t(r) = t. Analogously to pQ the random process Y(t), t > may be considered as a model 
of the anomalous diffusion on the interval [—1, 1] with reflecting edges. It is clear that if 
the curve T coincides with some straight line then the curve 7 is constructed in accordance 
with the law of geometrical optics: falling angle equals to reflecting angle. 



4 Analytical Representation of Reflected Diffusion 
Process and Its Properties 

Additionally to geometrical representation of the process Y(t), t > consider its analytical 
representation by the functions /:/?—>• [—2, 2), g : [—2, 2) — > [—1, 1], 

f(u) = (u + 2) /mod 4 -2, 

{u, —1<u<1, 
2 — it, 1 < u < 2, 
— 2 — u, -2<u< -1, 

Y(t)=g(f(y(t))). (4.1) 

Here u/mod A = A{u/A}, A > 0; {z} is the fractional part of a real number z. It is known 
from group theory that for u, v G R, A > 

(u/mod A) /mod A = u/mod A, {u + t>) /mod A = [u/mod A + v /mod A) /mod A. (4.2) 
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Define on the half-interval [—2, 2) binary operation " © " and unary operation of an inverse 

" e " : 

u®v = f(u + v), eu = /(-«), u,ve[-2, 2). 

Prove that these operations on the half-interval [—2, 2) form commutative group C with the 
summation © and the inverse and the unit : 

u®v = v®u, m©(© m ) = 0, 0©t>) ®w = u® (v@w), U,V,WEC. 

Really using ()4.2|) obtain 

u © v = f(u + v) = (u + v + 2) I mod 4-2 = (w + m + 2) /mod 4 - 2 = f(v + u) = v © u, 

u © (eu) = (u+ (Qu) + 2) /mod 4 - 2 = (u + (-« + 2) /mod 4-2 + 2) /mod 4 - 2 = 

= (u-u + 2) /mod 4 - 2 = /(0) = 0, 

(u@v)@w) = [f(u + v) + w + 2}/mod 4-2 = [(m + v + 2) /mod 4 - 2 + w + 2] /mod 4 - 2 = 

= (u + v + w + 2) /mod 4 - 2 = [u + (t> + w + 2) /mod 4 - 2 + 2}/ mod 4 - 2 = 

= [u + f[v + w) + 2]/mod 4 - 2 = m © (v © w). 

Prove that / : (i?, +) — > (C, ©) is a homomorphism of the additive group i? of real 
numbers onto the group C : 

/(« + «) = /(«)e/(«), f(-u) = ef(u). 

Using (|4.2|) obtain 

/(«) © f{v) = f(f(u) + f(v)) = {(u + 2) /mod 4 - 2 + (v + 2) /mod 4 - 2 + 2]/mod 4 - 2 = 

= [u + v + 2) /mod 4 - 2 = /(u + u), 
the equality f(—u) = Qf(u) is a corollary of the inverse definition. 

5 Reflection Formula for Diffusion Process 

Denote by pt = Pt(u), itt = ^t{u), Pt = Pt{u) distribution densities of random variables (r.v.) 
y(t), Y(t), f(y(t)) accordingly. Using the formula (|4.1|) and the graphic of the function / 
find: 

oo 

P t (u)= ^Pt(u-4k), ue[-2,2); P t (u) = Q, u#[-2,2). (5.1) 

v. f(v)=u k=—oo 

It is easy to prove from jHl chapt. 17, §6, lemma 1] that the row (|5.1jl converges. 
Define auxiliary function 

oo 

Pt(u) = / p t (u — Ak), — oo < u < oo. 
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The function Pt{u) coincides with Pt{u) for u G [—2, 2) and has the period 4 and is sym- 
metric. So from the formula (|4.1|) and the graphic of the function g obtain 

oo 

n t (u) = P t (u) + P t (u + 2)= ^2 Pt (u-2k), ue[-l, 1], ir t (u) = 0, u g [-1, 1]. (5.2) 

fc=— oo 

Suppose that r.v. z is distributed on [—1, 1], its density fi(v) has continuous derivative, 

^ = 0, , = ±1 
dv 

and r.v. z doesn't depend on random process y(t) then the formula (J5.1j) allows to calculate 
distribution density n^(u) of the random process Z(t) = g(f(y(t) — 1 + z)), t > in the 
moment t > : 

7r t 2 ( M ) = y (F t (u-u)+P t (u + 2 + u))/i(T;)dT;. (5.3) 

Remark that the formulas ()5.2j) . (|5.3j) which give distribution of diffusion reflected process 
are analogous to reflection method formulas which give solution of wave equation for finite 
string with fixed edges [U chapt. Ill, §13, points 5, 6]. 



6 Convergence of Reflected Diffusion Process 
Distribution to Uniform Distribution 

As / is the homomorphism of the additive group R onto the group C, then for < t\ < t 2 

f(y(t2)) = f(y(ti))®f(y(t 2 )-y(t l )). 

So random process f(y(t)), t > with independent and stable distributed increments is 
homogenous Markov process with state set [—2, 2) and symmetric distribution density of a 
transition from the state u into the state v during the time t 

q t (u, v) = P t (v -u), -2 < u, v < 2. 

Denote Q t = mf{q t (u,v), —2<u, v < 2} then from (j5.1j) obtain 

< Q t < i t > 0. (6.1) 



Suppose that 



^>h^J; 2 2) 2) ' (-) 



71 = 7t(u) 



1/2, u E [-1, 1], 
0, u £ [-1, 1]. 



For function cp on (— oo, oo) define the norm ||<^|| = sup{|<^(u)|, — oo < u < oo} and denote 
[z] the integer part of a real number z. 
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Lemma 1. For arbitrary h > 0, t > h : 

lk-7r|| <2(l-4Q h ) k - 1 \\P h -P\\ , k=[t/h}. 
Proof. As transition density qt(u,v) is symmetric then the formula ()6.2|) gives 
f 2 1 f 2 

/ P(v)q t (v,u)dv = - / q t (u,v)dv = P(u), —2<u<2. 

J -2 ^ y_2 

If Ag 4 = g t — > then the formulas ()6.1jl . (|6.4|) give 



(6.3) 



(6.4) 



|P m - P|| = sup 



Pi{v)q t (v,u)dv - / P(v)q t (v,u)dv 



-2<u<2 



sup 



2 z-2 

Pi(f) Ag t (v, — / P(v)Aq t (v,u)dv 

2 J-2 



-2 < u < 2 



and for t > 



|P t+ i-P|| < ||Pi - P||sup (J Aq t (v,u)dv, -2<u<2^j = \ \P X - P\ |(1 - AQ t ). (6.5) 



The formulas ()6.4|) . ()6.5|) allow to obtain by induction 

lin--P|| < (i-4g 1 ) fc - 1 ||p 1 -p|| , fc>i. 



(6.6) 



Using the formulas (j6.1|) . ()6.5j) it is easy to generalize the inequality ()6.6|) onto arbitrary 
t > 1 

HP-PH < (l-4g 1 ) fc - 1 ||P 1 -P|| , k=[t\. (6.7) 
The formulas ()5.2|) . ()6.7|) give 

||vr t -7r|| <2||P-P|| <2(l-4gi)^ 1 ||Pi-P|| , fc = [t], t>l. (6.8) 

The inequality ()6.8j) may be rewritten for arbitrary positive h in the form (jfi.HJl . □ 



7 Diffusion on Interval [— r, r] 

Suppose that m > 0, r = m _1//a and consider Markov process Y r (t), t > : 

Y r (t)=rg(f(y(t)/r)). 

The process V r (t), t > is obtained from the process ?/(£), t > by reflections from edges of 
the interval [— r, r]. Denote by 7r fi r = 7r t r (-u) distribution density of r.v. It is obvious 

that 7r t) ! = n t . Introduce normalized r.v. 

Wt , r = ^- = g(f(y(t)/r)),t>0 



with distribution density ip t ,r(u) = rTT tj , r (ru). 



Lemma 2. For K = [tin], tm > 1 



K, r(«) - ir(u/r)/r\\ < 



2(1-4Q 1 ) R - 1 \\P 1 -P 



(7.1) 



Proof. The definition of stable distribution (see chapt. 17, §5]) gives that for t > r.v.'s 
y(tm), y{t)/r coincide by distribution. So r.v.'s W t>r , g(f(y(t m ))) coincide by distribution 
too and from (J6.8)) 



7T 



\\A,r - 7r|| = \ Wm 

The formula (|7.1|) is proved. 

The formula (j7.1J) may be rewritten conditionally as 

r 1 (a,r,l) = r 1 (a,l,l)r a 



< 2(1 -4Q 1 ) K - i ||P 1 -P|| , K=[tm], tm>\. 



□ 



(7.2) 



8 Numerical Experiment 

Here analytical results of last paragraphs are compared with results of a numerical experi- 
ment. In this experiment a closeness of probability densities ir t , n is investigated for different 
t > 0. For this aim we imitate independent r.v.'s which coincide with r.v. Y(t) and so with 
r.v. 5f(/(t 1 / a y(l))) by distribution. R.v. y{\) is imitated approximately by normalized sum 



V i + . . . + V N 



of independent identically distributed r.v. v\ 



,...,v N , 



t~ a \t\~ a 

P{vi>t) = —, P{v x < -t) = £>L 

Using M independent realizations of r.v. g(f(t l l a y{l))) it is possible to calculate frequences 
Sj(t), j = 0, . . . , 9 of these realizations scorings into the sets 



10' 10 ' J 



2? 

1 + — , 1 
10' 



j =9. 



We calculate quantities 



S(t) = lOAfj^ ( s At) 



which characterize (analogously to x 2 -statistics) deviations of Y(t) distribution densities 
from uniform density for different t > 0. 
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t 

V 


S(t) for o— 1 9 

kJ \V J IV .71 KJj J. . U 


,9 (71 for (7—1 95 


S(t) for o— 1 99 


o ni 


91 69 31 


7867 49 


6584 26 


o D2 


4055 39 


3053 09 


2537 75 


n 03 


2200 99 


1 298 1 5 


1017 17 

lull .11 


n D4 


1 061 1 3 


609 86 


347 35 


0,05 


488.14 


289.15 


185.43 


0,06 


220.69 


177.32 


122.32 


0,07 


135.03 


52.97 


39.26 


0,08 


102.79 


39.37 


13.26 


0,09 


54.89 


14.55 


12.78 



Table 1 



Results of these calculations represented in the table 1 show how fast do distributions of Y(t) 
converge to uniform if the time t increases. Qualitative coincidence of numerical experiment 
results with the formula ()6.3|) is demonstrated. 



9 Diffusion with Periodical Initial Conditions 

Diffusion process with periodical initial conditions origines for example in fuel mixing at 
straight flow engine [2]. To model this process take natural n and define markov process 

Z n {t)=g{f(y{t)-l + z n )), t>0 

where r.v. z n has uniform distribution on finite set 

2k + I 

= : k — 0, 1, . . . , n — 1 

n 

and z n , y(t) are independent. Then random process Z n (t), t > may be considered as 
anomalous diffusion on interval [—1, 1] but with periodical initial conditions 

P(Z n (0) = -1 + 8) = -, SEl n . (9.1) 

n 

Denote by IT, „ = IT, n (u), Pt, n = Pt, n(u) distribution densities of r.v.'s Z n (t), f(y(t) — 
l + z n ), t>0. 

Lemma 3. The following formula is true 



TT / X 1 ^ ( 2k+\\ 

n t,»M Hn [U+1 — " J • 

Proof. Analogously to the formulas ()5.1j) - ()5.3|) obtain 

1 1 °° 

Pt, n{u) = - }2 p t( u + 1 - s ) = ~ J2 J! Pt(u-4k + l-s), -2<u<2, 



(9.2) 



n * — ' n 

sGl n sdl n k=—oo 



7 



n t ,„0)= Pt,n(v) = -J2(P t (u+l-s) + P t (u + l + S )), ~1<U<1 

v. g(v)=u sGl n 

then 

n t , n (u) = - V jh ( u - — ) , -1<u<1. (9.3) 

k=-oo v 7 

The formula ()9.3|) leads to 

n t> „(«) = n t , „ (u + , -i<w<i-i (9.4) 

Calculate now the function 

2 

n t , n (it), -l < u < -l + -. 

n 

For this aim take 

1 1 1 

u = w — 1 H — , < u; < — . 

n n n 

Then in an accordance with (|5.2|) . (|9.3|) 

n 4 , n H = > — < tW < -. 9.5 

n n n 

The formulas flO}, (jHSJ) lead to the equality (Q. □ 

The equality ()9.2j) means that the diffusion (normal or anomalous) on the set [—1, 1] 
with periodical initial conditions (J9.1|) and reflecting edges leads to the same result as a 
diffusion on isolated (by reflecting edges) subsets 



2k + 1 1 2k + 3 1 
-1 + , -1 + + - 

n n n n 



k — 0, . . . , n — 1 



of the set [—1, 1]. 

Remark that the equality is true for each process y(t) with independent and symmetrically 
distributed increments. 

Theorem 1. For tn a > 1, L = [tn a ]: 

||n t , n -7r|| <2(1-4Q 1 ) L - 1 ||P 1 -P|| . (9.6) 

Proof. The statement of the theorem is obtained directly from the equality ()9.2|) and the 
formula ()7. 1|) in which r equals to 1/n. □ 

The inequality (j9.6J) may be interpreted as a decreasing of characteristic time into n a 
times in model of anomalous diffusion with n— periodical initial conditions ()9.1|) : 

T n (a, 1,1) = T 1 (a,l,l)/n°. 

This result may be easily spread onto general case when r.v. Z n (0) has distribution density 
with continuous derivative r n {u) which satisfies periodicity conditions: 

r n (u) =r n (u+-) , -1 < u < 1 



n I n 



and symmetry conditions 



r n (-l + - - v) = r n ( -1 + - + v ) , < v < - 
n \ n J n 



0, v = ±l. 



and boundary condition 

dr n (v) 
dv 

It is a generalization of [2. results from normal onto anomalous diffusion. 

10 Mult i- Dimensional Diffusion with Periodic Initial 
Conditions 

Multi-Dimensional Anomalous Diffusion. Results of last paragraph may be generalized 
onto two-dimensional case if the interval [—1, 1] is replaced by square [—1, 1] x [—1, 1] and 
process of anomalous diffusion on square with reflecting boundaries is considered as two 
independent processes of one dimension anomalous diffusion on the intervals [—1, 1] with 
reflecting edges. Disintegrate the unit square by rectangular network with n 2 equal squares. 
Suppose that initial state of two-dimension diffusion process on unit square has uniform 
distribution on the set of these n 2 squares centers. Then all one-dimension estimates of 
convergence rate are spread onto two-dimension case. 

Denote all one-dimension distributions from last paragraph by indexes characterized num- 
bers of appropriate coordinates j = 1,2: 

P (j) p? = Pi ] {uj), tt<0 = n«l = n«l( Uj ), j = 1,2. 

In accordance with the model of two-dimension diffusion (as a pair of independent one- 
dimension diffusion models) obtain: 

7T = 7r(«i, U 2 ) = 7r (1) («i)7r (2) (« 2 ), Rt,n = ^-t, n(«l, U 2 ) = Tl^U^Tlf^U^ . 

From the inequality flSEJ for L = [tn 1 / a ], tn 1 ^ > 1, j = 1, 2 

||ITf n - vr^H < 2(1 - 4Q 1 ) L ~ 1 \\P^ - P<fl\\ = , . (10.1) 

Define the norm ||$|| = sup{|$(wi, u 2 )\, — oo < u%, u 2 < oo} of the function $ on a plane 
and put A t , n = A^ n = A^ n . 

Theorem 2. For L = [tn 1/a ], tn 1/a > 1 

||n t , n -7r|| < A t>n (l+A t>n ). (10.2) 

Proof. From triangle inequality obtain 

I in, „-7r|| = | |n« nf n - ^ V 2 ) \ \ < \ |<> nf n - ^{\\ + 1 \^n® - \ \ < 
< lin^H ||n« - + li^ll ||nf„, - ^\\ < (\\^\\ + - nf n ||) ||<> - 

lina-««|| ^ g + IN«-na,|) ||nft--«|| + i ||nR-<«|| 

Using the formula (|10.1|) obtain the inequality (|10.2|) . □ 
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As Af t „ — > 0, t — > oo so convergence rate of two-dimension case (J1U.2)) is analogous 
to one- dimension case ([9.6)1 . This result may be easily spread onto multi-dimension case. 
Specifics of this multi-dimension diffusion model is that as one-dimension diffusion compo- 
nents so their initial conditions are independent. 

Multi-Demensional Normal Diffusion. Consider /c-dimension normal diffusion model 
in which U t n (yi, ■ ■ ■ , Uk) is distribution density in the moment t at the point (yi, . . . , y^) G 
[-1,1]*: 



k Q2 



d Q 2 

d t -J2dy]J U tyn (y u ...,y k )=0, 
An t>n (±l,...,±l,y i ,±l,...,±l)=0, -l< Vj <l, j = l,...,k, 

_^ m k 

n , n (yi, ■ ■ ■ , Vk) = git + a ( Jl > II cosTmj r y r > 0, m < oo. 

ji,-Jfe=i r=1 



Then 



TL t>n {yi,---,yk) = ^ + ^ JJcosvrnj r ?/ r exp -vrVt^A;| 



iiv-Jfc=l r=l \ i=l 



|a(l,...,l)| >0 



||II tin (yi, . . . ,2/ fc ) - ^|| ~ |a(l, . . . , l)|exp(-7r 2 fcn 2 t), i -> oo. 

So it is possible to choose n-periodic (by each coordinate) initial conditions that 

T n (2,l,k)=T 1 (2,l,l)/kn 2 . 
The paper is supported by RFBR, projects 03-01-00512, 03-07-07-90334v. 
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